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ABSTRACT 

WASP- 14b belongs to a class of highly UTadiated hot Jupiters (Teq = 1866K) with a mass of 7.3 ± 0.5 My and 
a radius of 1.28 ± 0.08 Rj. With a mean density of 4.6 gcm'^, this transiting planet is one of the densest known 
to date (Joshi et al. 2009), for planets with periods less than 3 days. We present analytic light-curve models 
for three Spitzer secondary eclipses, Keplerian orbital model, estimates of infrared brightness temperatures, 
and constraints on atmospheric composition and thermal structure. Although extremely irradiated, WASP- 14b 
does not show any distinctive evidence of a thermal inversion, based on existing observations. In addition, the 
observations indicate low day-night redistribution { fr < 0.25) and are consistent with chemical compositions 
at nearly solar abundances. We confirm a significant eccentricity of e = 0.087 ± 0.002 and refine other orbital 
parameters. 

Subject headings: eclipses — planetary systems — planets and satellites: atmospheres — planets and satellites: 
individual (WASP- 14b) — techniques: photometric 



1. INTRODUCTION 

The Spitzer Space Telescope (Werner et al. 2004) is the 
most widely used facility for characterizing extra solar plan- 
ets. Out of space telescopes, Spitzer systematics are among 
less complex and are very well-studied and modeled, provid- 
ing an invaluable resource for accurate exoplanet characteri- 
zation (Seager & Deming 2010). This has enabled the mea- 
surement of tens of atmospheres, using the detection of pri- 
mary and secondary eclipses as the most prolific method of 
investigation to date. 

The planet-to-star flux ratio is enhanced in the infrared due 
to the planet thermal emission and simultaneous drop in the 
stellar emission, enabling high-precision photometric mea- 
surements. Combining several secondary-eclipse observa- 
tions measured in broad Spitzer bandpasses with the Infrared 
Array Camera (IRAC, Fazio et al. 2004), a low-resolution 
dayside spectrum from the planet can be reconstructed, re- 
vealing key atmospheric and physical parameters. These mea- 
surements can further be used to constrain atmospheric com- 
position, thermal structure, and ultimately the formation and 
evolution of the observed planet. WASP- 14b represents an in- 
triguing object for such an analysis, having characteristics not 
so common for close-in, highly irradiated giant planets. 

WASP- 14b (Joshi et al. 2009) was discovered as a part of 
the SuperWASP survey (Wide- Angle Search for Planets, Pol- 
lacco et al. 2006; ColUer Cameron et al. 2006, 2007). Photo- 
metric and radial velocity observations revealed a planetary 
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mass of 7.3 ± 0.5 Mj and a radius of 1.28 ± 0.08 Rj. Its 
density (p = 4.6 gem ■') is significantly high compared to typ- 
ical hot Jupiter densities of 0.34-1.34 gem""* (Loeillet et al. 
2008a). The planet is also very close to its star (semi-major 
axis 0.036 ± 0.001 AU), and has a significant orbital eccen- 
tricity of 0.087 ± 0.002 (this work). 

Detailed spectroscopic analyses of the stellar atmosphere 
determined that the star belongs to the F5 main-sequence 
spectral type with a temperature of 6475 ± 100 K and high 
lithium abundance of log N(Li) = 2.84 ± 0.05. F-type stars 
with this temperature should have depleted Li, being close to 
the Li gap or 'Boesgaard gap' (Boesgaard & Tripicco 1986; 
Balachandran 1995). However, the high amount of Li and a 
relatively high rotational speed of vsin(/) = 4.9 ±1.0 kms"' 
indicate that WASP- 14 is a young star. Comparing these re- 
sults with models by Fortney et al. (2007) for the range of 
planetary masses and radii led Joshi et al. (2009) to constrain 
the age of the system to 0.5-1.0 Gyr 

Joshi et al. (2009) also discuss the high eccentricity of the 
planet. Because WASP- 14b has a very small orbital distance, 
probable scenarios for such a significant eccentricity would 
be either that system age is comparable to the tidal circular- 
ization time scale or there is a perturbing body. 

Husnoo et al. (2011) performed long-term radial velocity 
measurements with a goal of discovering or rejecting the pres- 
ence of a third body. They refine the orbital eccentricity value 
(e = 0.088 ± 0.003) and spin-orbit angle. They argue that this 
planet has undergone some degree of orbital evolution, but it 
is still subject to strong tidal forces. This establishes a new 
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lower limit for the semimajor axis at which orbital eccentric- 
ity can survive tidal evolution for the age of the system. They 
state that since there is no observable unambiguous trend in 
residuals with time, there is no firm evidence for an existing 
planetary companion. 

Using the Spitzer Space Telescope, we observed three sec- 
ondary eclipses at 3.6 jam, 4.5 jam, and 8.0 ^.m. We present 
analytic light-curve models that incorporate corrections for 
systematic effects with a new approach to removing them with 
the Stevenson et al. (2011) pixel sensitivity mapping tech- 
nique, a Keplerian orbital model, estimates of infrared bright- 
ness temperatures, and constraints on atmospheric composi- 
tion and thermal structure. 

In Section 2 we describe our observations. Section 3 dis- 
cusses data reduction procedures. Sections 4 presents our 
photometry and Section 5 discusses the modeling techniques 
and results of each dataset. Section 6 presents constraints on 
the orbit of WASP- 14b, and Section 7 reveals the atmospheric 
structure and composition. In Section 8 we discuss our results 
and in Section 9 we present our conclusions. 

2. OBSERVATIONS 

The Spitzer IRAC instrument observed two events; one at 
3.6 i^m in 2010 March (Knutson's program 60021) and one si- 
multaneous observation in two wavelength bands (4.5 and 8.0 
(im) in 2009 March (Harrington's program 50517). Observa- 
tion at 3.6 \i.m (channel 1) was made in subarray mode with 2 
s exposure, while observations at 4.5 and 8.0 |i.m (channels 2 
and 4) were made in stellar mode (2x2/12) with pairs of 2 s 
frames taken in the 4.5 ^m band for each 12 s frame in the 8.0 
|im band. This mode was used to avoid saturation in channel 
2. 

TABLE 1 

Observation Information 



Channel 


Observation 


Start Time Duration Exposure Number of 




Date 


[JD] [s] 


Time [s] 


Frames 


Main Science Observation 


Chi 


2010-03-18 


2455274.4707 28055.4 


2 


215 


Ch2 


2009-03-18 


2454908.8139 19998.7 


2x2 


2982 


Ch4 


2009-03-18 


2454908.8139 19998.7 


12 


1481 


Pre-observation 


Ch2+4 


2009-03-18 


2454908.7877 2019 


2 


213 


Post-observation 


Ch244 


2009-03-18 


2454909.0455 367 


2x2,12 


10 



We have pre- and post-calibration frames for the 4.5 and 8.0 
jxm observation. Prior to the main observation, we exposed 
the array to a relatively bright source. That quickly caused the 
saturation of the detector, and produced less-steeper increase 
during the main observation (Section 5.2). During the short 
post-eclipse observation, frames of the blank sky allowed us 
to check and demonstrate that there are no warm pixels in the 
aperture. The Spitzer pipeline version used for 3.6 jxm ob- 
servation is S. 18. 14.0 and for 4.8 and 8.0 |xm observation is 
SI 8.7.0. The start date of each observation, duration, expo- 
sure time and total number of frames are given in Table 1 . 

3. DATA REDUCTION 

3.1. Background 

Our analysis performs Photometry for Orbits, Eclipses and 
Transits (POET). It produced light curves from Spitzer Basic 
Calibrated Data Frames (BCD), fits models to the light curves, 



and assesses uncertainties. The derived parameters constrain 
a separate orbital and atmospheric model. In this section we 
give a general overview of POET's routines. The following 
sections will explain in detail those relevant for our analysis. 

Each analysis starts with identifying and flagging bad pix- 
els in addition to the ones determined by the Spitzer bad pixel 
mask (see Section 4). Then we perform centering. Due to the 
^0.1% relative flux level of secondary-eclipse observations 
and Spitzer' s relative photometric accuracy of 2% (Fazio et al. 
2004), we apply a variety of centering routines, looking for 
the most accurate one. Tliree methods are used to determine 
the PSF center precisely: center of light, two-dimensional 
Gaussian fitting, and least asymmetry (see Lust et al. 2012). 
We then apply 5 x -interpolated aperture photometry (Harring- 
ton et al. 2007), where each image is re-sampled using biUn- 
ear interpolation. This allows the inclusion of partial pixels, 
thus reducing pixeUation noise (Stevenson et al. 2011). We 
subtract the background within a specified annulus centered 
on the star, and during testing of the large range of aperture 
sizes, we discard the frames with bad pixels in the photometry 
aperture. 

To explore the parameter space and estimate uncertainties, 
we use our Markov Chain Monte Carlo routine (MCMC, see 
Section 3.3), where we model all known systematics and 
the eclipse event simultaneously, using at least lO'' iterations 
(Stevenson et al. 2010; Campo et al. 2011). To test adequate 
sampling, we ran four independent MCMC chains. To ensure 
the convergence of all chains, we run as many iterations as 
necessary until the Gelman & Rubin (1992) test for all free 
parameters drop below 1%. Our MCMC routine can model 
events separately or simultaneously, sharing parameters such 
as the eclipse midpoint, ingress/egress times or duration. 

To find the best model, we test all systematic models known 
in the literature and some of our own. We use our newly de- 
veloped (Stevenson et al. 201 1) BiLinearly-Interpolated Sub- 
pixel Sensitivity (BLISS) technique to model intrapixel sensi- 
tivity variation at the sub-pixel level (see Section 3.2). Using 
both bilinear (BLI) and nearest-neighbor (NNI) interpolation 
and over 1000 knots (the control points), we map the pixel 
surface. This method does not require any free parameters 
and quickly converges. 

To determine the best aperture size, we choose the solution 
with the smallest standard deviation of normalized residuals 
(SDNR). The best model for a given aperture size is deter- 
mined by applying the Bayesian (BIC) and Akaike (AIC) in- 
formation criteria, which compare models with different num- 
bers of free parameters (see Section 3.3) . 

We also base our decision upon the level of correlation in 
the residuals (Campo et al. 2011; Nymeyeret al. 2011). Less- 
correlated models are better, if there is no significant differ- 
ence in AIC and BIC (see Sections 3.3 and 5.1). 

Finally, we report Julian Date mid-times in both BIDy-j^^ 
(Coordinated Universal Time) and BJD^^ (BlDj^g, Barycen- 
trical Dynamical Time), calculated using the Jet Propulsion 
Laboratory (JPL) Horizons system, to facilitate handling dis- 
continuities due to leap seconds and to allow other researchers 
to compare mid-times of the eclipses easily. 

3.2. Modeling Systematics 

Modeling systematics is critical to accurate atmospheric 
characterization. So far, there are two major discrepancies 
presented in the literature. Tinetti et al. (2007) find a deeper 
transit of HD189733b flian Desert et al. (2009), which in- 
dicates a larger optical depth and the presence of water va- 
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per higher in the atmosphere. Knutson et al. (2009b), re- 
analyzing 8.0 \xm observations of HD149026b, by applying 
a different set of systematic models, find a range of eclipse 
depths from 0.05 - 0.09%, preferring the lower values as their 
best. This value is half that of Harrington et al. (2007) (see 
Stevenson et al. (2011), for our latest analyses of this data 
set). To avoid similar and future errors, our team is heavily 
focused on modeling systematics, testing all of the known an- 
alytic model components. 

Spitzer's IRAC has four broad photometric bandpasses: 3 .6, 
4.5, 5.8, and 8.0 \im (channels 1, 2, 3 and 4). All four chan- 
nels exhibit large systematic errors on the scale of our re- 
quirements. The 3.6 and 4.5 |i.m bands use InSb detectors, 
and the 5.8 and 8.0 \im bands Si:As detectors. InSb detec- 
tors are known to show variable response in quantum effi- 
ciency depending on the centroid location within the pixel. 
This intrapixel sensitivity is greatest at pixel center and de- 
clines toward the edges of the pixel, typically by several per- 
cent (Morales-Calderon et al. 2006). It is also not symmetric 
about the center of the pixel and the amplitude of the effect 
varies from pixel to pixel. Over the total duration of the obser- 
vation, the position varies by several tenths of the pixel width 
and since the motion of the stellar center goes through this re- 
gion frequently, this systematic is adequately sampled during 
a single eclipse observation. One minimizes this effect by ob- 
serving with fixed pointing (Reach et al. 2005; Charbonneau 
et al. 2005; Harrington et al. 2007; Stevenson et al. 2010). 

The other systematic, the most prominent and well studied, 
is the time varying sensitivity, called the ramp effect. This 
is a baseline effect which may affect all the channels, but is 
most pronounced at 8.0 [im (Charbonneau et al. 2005; Har- 
rington et al. 2007). It is observed as apparent increase in 
flux with time, which is believed to be caused by charge trap- 
ping. An eclipse is superimposed on this baseline and is eas- 
ily separated, but not without adding uncertainty to the eclipse 
depth. After many hours of observations this effect reaches a 
plateau, hence a way to stabilize it is to point the array to a 
bright source prior to the main observation. This causes satu- 
ration and consequently produce a flatter ramp (Knutson et al. 
2009b; Seager & Deming 2009). The systematic is particu- 
larly important for weak eclipses, where an incorrect choice 
of ramp can produce an incorrect eclipse depth. 

Although each type of systematic, due to different type of 
alloy on the detectors, is strong in different set of channels 
(ramp effect channels 3 and 4, intrapixel effect, channels 1 
and 2), it was recently reported that both systematics were 
seen in other channels (Stevenson et al. 2010). That is why 
we model all variations in all channels. 

To model ramp effect, we use different forms of exponen- 
tial, logarithmic, and polynomial models used in the literature, 
and some of our own, from more than 10 different equations 
implemented in our pipeUne (see the full set of equations in 
Stevenson et al. 2011). For each channel we test all ramp 
models, assessing the best one by using criteria listed at the 
end of Section 3.3. 

Intrapixel variability, which may vary up to 3.5% within the 
pixel, is modeled with our new intrapixel technique, caUed 
BLISS mapping (Stevenson et al. 2011). This technique, op- 
posed to the widely used two-dimensional polynomial fitting 
(Knutson et al. 2008; Machalek et al. 2009; Fressin et al. 
2010) and the mapping technique developed by Ballard et al. 
(2010), maps the spatial variability on a subpixel- scale grid 
using bihnear interpolation, and shows significant improve- 
ment in both SDNR and BIC compared to other methods. See 



Stevenson et al. (201 1) for comparison tests. 

This method first divides the observed flux by eclipse and 
ramp models, assuming that any residual fluxes are related to 
the pixel position. Then it creates a rectangular xy grid of 
knots at subpixel resolution, which are used as control points 
in the interpolation. The flux value of a particular knot is the 
mean of the residual fluxes assigned to that knot. 

We use over 1000 knots to map the pixel surface and two in- 
terpolation calculations, bilinear (BLl) and nearest-neighbor 
(NNI), to generate the sensitivity map (see Section 5.1). 
Nearest-neighbor assigns each point the value of its nearest 
knot, and does not consider the values of neighboring points, 
producing a piecewise-constant interpolant. See Stevenson 
et al. (201 1) for more detail about this technique. 

The crucial part in this analysis is determining the knot 
spacing (bin size), which gives resolution in the position 
space. The bin size must be smaU enough to catch any smaU- 
scale variation, but also large enough to ensure no correla- 
tion with the eclipse fit (see Section 5.1). Assuming accurate 
centering, BLl should always outperform NNI. The bin size 
where NNI outperforms BLl is indication of the centering pre- 
cision, and determines the bin size for the particular data set. 
If NNI always outperforms BLl, that is the indication that in- 
trapixel variabiUty is very weak and intrapixel modeling is 
unnecessary. 

3.3. Modeling Light Curves and the Best Fit Criteria 

To find the best model, for each aperture size we systemati- 
cally explore every combination of ramp model and intrapixel 
sensitivity model. The final light curve model is expressed as, 

F{x,y,t)=FMt)M{x,y)E{t), (1) 

where Fix.y.t) is the aperture photometry flux, F^ is the 
constant system flux outside of the eclipse, R{t) is time- 
dependent ramp model, M{x,y) is the position-dependent in- 
trapixel model and E{t) is the eclipse model (Mandel & Agol 
2002). 

To estimate uncertainties we use our MCMC routine with 
the Metropohs-Hastings random walk algorithm with at least 
10* iterations, to ensure accuracy of the result. This rou- 
tine simultaneously applies eclipse parameters and Spitzer 
systematics, exploring the parameter phase space, determin- 
ing uncertainties and correlations between the parameters. It 
leaves depth, duration, midpoint, system flux, and ramp pa- 
rameters free. Additionally, the routine can model multiple 
events at once, sharing the eclipse duration, midpoint and 
ingress/egress times. These joint fits are particularly appro- 
priate for channels observed together (see Campo et al. 2011 
for more details about our MCMC routine). 

To avoid fixing any model parameter during MCMC, we 
use the Bayesian prior (Gelman 2002). This is particularly 
relevant for noisy or low S/N datasets where some parameters 
are not well constrained by the observations, as is the case for 
ingress and egress times. For them we use informative priors, 
taken from the literature. 

Photometric uncertainties used in our analyses are derived 
by fitting an initial model with a minimizer and re-scaUng 
it to reduced of 1 , since Spitzer pipeline uncertainties are 
generally over-estimated (Harrington et al. 2007). This scal- 
ing factor is proportional to SDNR of the models. Upon 
rescaUng, we re-run the least-squared minimizer with the new 
uncertainties, and then MCMC. If MCMC finds a lower 
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Fig. 1. — Raw (left), binned (center) and systematics-corrected (right) primary transit liglit curves of WASP-14b at 3.6, 4.5 and 8.0 |j.m. Tlie results are 
normalized to the system flux and shifted vertically for the purpose of comparison. The colored Hues are best-fit models, the black curves omit their transit 
models, and the error bars are la uncertainties. As seen in the binned plots of wa014bs21, wa014bs41, a linear function accurately fits the time-dependent 
systematic without the need to clip any initial data points. wa014bsl 1 fits better if initial data points are clipped (see Section 5.1). 



value than the minimizer, we re-run the minimizer starting 
from the MCMC's best value. The minimizer will find an 
even better x^- We then restart the MCMC from the new min- 
imizer solution, and run 10*' iterations to derive uncertainties 
for each of the parameters (Campo et al. 201 1) 

Our measures of goodness of fit are SDNR, BIC and AIC 
values (Liddle 2007): 



BIC = x +'tlnA^, 



AlC = x +2k, 



(2) 



(3) 



where k is the number of free parameters, N is the number of 
data points. These criteria penalize additional free parameters 
in the system, with better fits having lower values. To appro- 
priately compare BIC and AIC values for the same aperture 
size, and determine the best fit, we use the same uncertainties 
for each dataset, and model all combinations of ramp models 
and intrapixel model. SDNR values are used to compare dif- 
ferent aperture sizes. The lowest value defines the best aper- 
ture size. 

Equally important is the correlation in the residuals (see 
Section 5.1). We plot and compare the theoretical 1 / \/N scal- 
ing with the RMS model residuals vs. bin size (Pont et al. 
2006; Winn et al. 2008). If there is no major deviation in 
those two curves, there is a lack of significant time-dependent 
photometric noise correlation in the final models. 

After MCMC is finished, we study parameter histograms 
and pairwise correlations plots, making sure that the phase 
space minimum is global. 

4. WASP- 14b PHOTOMETRY 

For our analyses, we used Basic Calibrated Data (BCD) 
frames generated in the Spitzer IRAC pipeline (Fazio et al. 
2004). The pipeline version used for each observation is given 
in Section 2. Our data reduction procedure started with apply- 
ing Spitzer s bad pixel masks and with our procedure for flag- 
ging additional bad pixels (Harrington et al. 2007). In each 



group of 64 frames and at each pixel position, we applied two- 
iteration outlier rejection, which calculated the frame median 
and the standard deviation from the median, and flagged pix- 
els that deviates more than 4a. Then we found the stellar cen- 
troid for the photometry, by using a 2D Gaussian fit to data in 
an aperture radius of four pixels. 

Light curves were extracted using 5 x -interpolated aperture 
photometry (Harrington et al. 2007), for every aperture radius 
from 2.25 to 4.25 pixels, in 0.25-pixel steps. The background 
is estimated and subtracted before the photometry, in an annu- 
lus centered on the source with inner and outer radii defined 
in Section 5.4. 

To calculate the Barycentric Julian Date (BJD) of each ex- 
posure we used the mid-exposure time of each frame, based 
on the UTCS-OBS value in the FITS header and the frame 
number. We performed our barycentric light-time correction 
using our own code and the coordinates of the Spitzer space- 
craft from the Horizons ephemeris system of the Jet Propul- 
sion Laboratory. The times are corrected to BJDj^g to remove 
the effects of leap seconds and relativity. 

5. WASP-14b SECONDARY ECLIPSES 

Here, we discuss each observation in detail, particularly fo- 
cusing on channel 1, due to the demanding analysis of that 
data set. In Subsection 5.1 we give our control plots, as an 
example of how we verify that our results are indeed the best 
solution for the particular data set. We present each channel 
separately and then a joint fit. We give our best model fits, 
with the particular ramp and intrapixel model used for each 
data set, corresponding eclipse depths and brightness temper- 
atures and joint best-fit eclipse parameters. Figure 1 shows 
our best fit eclipse light curves. In Appendix we give param- 
eters for the WASP- 14 system, derived from this analysis and 
from the literature. 

5.1. Channel 1 - 3 .6 \xm 

The observation in channel 1 lasted 7.8 hours, during which 
a good amount of the system baseline was observed before 
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and after the secondary eclipse event. Images that are taken 
at the start of the observation have larger pointing offset, due 
to the telescope settling prior to the eclipse observation. Dis- 
carding starting data in that case improves the fit. Some initial 
data points (^36 minutes, 1 100 of 13760 points) were there- 
fore removed from the fitting procedure. In Figure 2 we com- 
pare SDNR values for models with different ramps with and 
without exclusion of the initial data points. Models with ini- 
tial data points removed produce better fits with lower values 
for SDNR. 
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Fig. 3. — The best-fit eclipse depths as a function of photometry aperture 
size for wa014bsll. All four best ramp models are plotted (see bellow). 
The red point indicates the best aperture size for that channel. The eclipse- 
depth uncertainties are the result of 1000 MCMC iterations. The trend shows 
insignificant dependence of eclipse depth on aperture size (less than Icr). 
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Fig. 2. — Top: Channel 1 comparison of different ramp models with all ob- 
servational points included (no-preclip). Plot shows SDNR v.v. aperture size. 
A lower value indicates a better model fit. Bottom: Channel 1 comparison 
of different ramp models with 1 100 points excluded (preclip). 

Starting from the aperture size of 2.25 and continuing with 
an increment of 0.25-pixel, we tested all of the ramp mod- 
els (linear, rising, exponential, sinusoidal, double exponen- 
tial, logarithmic, etc.). Corresponding equations are listed in 
Stevenson et al. (2011). To determine the best solution we 
consider our best-fit criteria (see Section 3.3) and study the 
correlation plots. Most of the models produced obvious bad 
fits, so minimizer and shorter MCMC runs eliminated them. 
The best aperture radius is 2.75 pixels (see Figure 2, bottom 
pannel). We tested the dependence of eclipse depth on aper- 
ture radius (Anderson et al. 2010). The trend in some events 
may indicate a slightly imperfect background removal (see 
Figure 3). The effect is less than Icr on the eclipse depth. 

Intrapixel variability is the strongest systematic in Spitzer 
channels 1 and 2. The sensitivity BLISS map is given in 
Figure 4. The correlation coefficients between eclipse depth 
and BLISS map knots are given in Figure 5. As we stated in 
the Section 3.2 the most important variable to consider with 
BLISS is the bin size. It defines the resolution in position 
space and gives insight into the average behavior The posi- 
tion precision for channel 1, measured as RMS of the position 
difference on consecutive frames, shows significant difference 
between the x and y axes (see Figure 6, upper panel). We con- 
sidered a range of bin sizes for both BLI and NNI around the 



calculated precision. The best bin size for this data set was 
determined when NNI outperformed BLI. It is 0.004 pixels 
for jc and 0.01 fory. 
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Fig. 4. — Top: BiLinearly-Interpolated Subpixel Sensitivity (BLISS) map 
of wa014bsll. Redder (bluer) colors indicate higher (lower) subpixel sensi- 
tivity. The horizontal black lines defines the lower pixel boundary. Bottom: 
Pointing histogram. Colors indicate the number of points in a given bin. 

We also tested two-dimensional polynomial intrapixel mod- 
els (Knutson et al. 2008; Stevenson et al. 2010; Campo et al. 
2011): 

V]p(x,y) = piy^+P2x^+P3xy+p4y+p5x+l, (4) 

where x and y are the centroid coordinates relative to the pixel 
center nearest the median position and pi - ps sae free param- 
eters. Both SDNR and BIC show significant improvements 
with the new BLISS models over the polynomial models (see 
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Fig. 5. — Con'elation coefficients between eclipse depth and computed 
BLISS map loiots for wa014bsll. The correlation regions (in red) indicate 
that it is necessary to compute the BLISS map at each step of an MCMC 
routine, so that we could correctly asses the uncertainty on the eclipse depth. 

Table 2 . The improvement in BIC is expected, since neither 
BLI nor NNI require any free parameters. The improvement 
in SDNR confirms that interpolation offers a better fit to the 
data (see Stevenson et al. 2011 for other tests that compare 
polynomial and BLISS intrapixel models). 

TABLE 2 

Comparison BLISS And Best Polynomial Model 



Ramp Model 



BLISS 



Polynomial-Quadratic 





SDNR 


BIC 


SDNR 


BIC 


Channel 1 


No ramp 


0.003313 


12350.0 


0.0033853 


12593.2 


Linear 


0.003311 


12342.3 


0.0033852 


12588.5 


Sinusoidal 


0.003316 


12342.2 


0.0033855 


12590.5 


Quadratic 


0.003310 


12351.5 


0.0033850 


12597.3 


Channel 2 


Linear 


0.004472 


2971.9 


0.0045272 


2992.2 



In this channel, we detect photometric noise correlations in 
our models, with consistent residual variation observed in all 
aperture sizes, with all ramp models, and at every bin size 
(see Figure 6, bottom panel and Figure 7, left panel). Al- 
though channels 1 and 2 (3.6 and 4.5 |^m) both exhibit strong 
intrapixel variations due to the same material present on the 
detector (InSb), typically channel 1 exhibits larger sensitivity 
than channel 2 even if the star is placed in the middle of the 
pixel (Knutson et al. 2011). However, although the star falls 
close to the edge of the pixel, this time variation is unusually 
high, and cannot be modeled by any of our position- or time- 
dependent systematic models. This might indicate that this 
young, WASP- 14 star is still active. 

In our work on dozens of eclipses, we have often found the 
same Spitzer channel to behave differently at different times, 
even on the same star Our systematics removal algorithms 
correct the worst effects, which are consistent, but there is of- 
ten still some significant baseline scatter or oscillation. While 
one might expect certain kinds of stars to be relatively stable, 
recall that Spitzer reaches a ^ 0.01% sensitivity in these ob- 
servations, and stellar oscillations of this scale and at these 
wavelengths are not well studied. So, it is not fully clear 
whether these effects come from the observatory or the star. 
Since scatter and oscillation often persist during an eclipse, 
and since a change in planetary signal of the magnitude seen 
would generally mean an unbelievably dramatic change in the 
planet, we feel justified in treating the scatter or oscillation 
phenomenologically. If there is significant pointing offset at 
the start of the observation due to the telescope settling, we 




1.0010 

1.0005 

1.0000 

D.9995 

D.9990 

D.9985 

D.9980 

D.9975 

D.9970 
D.OOlO 
D.0005 
D.OOOO 
0.0005 
0.0010 
0.0015 




..,0 0.45 0.50 0..... 

Orbital Phase 12.24-day period) 

Fig. 6. — Top: Projected BLISS map of wa014bsll integrated along the 
-V (right) and y (left) axes. This position dependence variability is effectively 
fitted with our new BLISS mapping technique. Bottom: Residuals for the 
channel 1 observations (lower panel) display correlated noise both in and out 
of the eclipse, probably due to stellar activity. 
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Fig. 7. — Con'elations of the residuals for the three secondary eclipses of 
WASP-14b. Black color represents the RMS residual flux vs. bin size. The 
red line shows the predicted standard error for Gaussian noise. Large excesses 
above the red line indicate coirelated noise at that bin size (see Section 5.1). 

remove it by clipping large excursions and we choose the 
ramp model that fits best under BIC and AIC. Our MCMC 
accounts for any correlation between eclipse and model pa- 
rameters, and our RMS vs. bin size analysis accounts for any 
residual time correlation of the residuals. 

We performed 10^ MCMC iterations to determine the pa- 
rameter posterior distributions and their uncertainties. We 
used the Bayesian informative prior for the values of ingress 
and egress times, obtained from Joshi et al. (2009) and David 
R. Anderson (personal communication, 2010). All other pa- 
rameters (eclipse midpoint, eclipse duration, eclipse depth, 
system flux, and ramp parameters) were left free. 

TABLE 3 
Channel 1 Ramp Models 



Ramp Model 


SDNR 


BIC 


Eclipse 








Depth [%] 


No ramp 


0.003313 


12350.0 


0.1837 ± 0.0071 


Linear 


0.003311 


12342.3 


0.1867 ±0.0072 


Sinusoidal 


0.003316 


12342.2 


0.1932 ± 0.0071 


Quadratic 


0.003310 


12351.5 


0.1896 ±0.0098 



Considering all the above criteria (see also section 3.3), we 
selected 4 ramp models that produce almost identical SDNR 
and BIC values. We concluded upon studying output results, 
histograms and correlation plots that there is no single best so- 
lution for this data set. Because quadratic and sinusoidal ramp 
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models generally produce deeper eclipse depths, for the final 
atmospheric analyses we selected two extreme final results 
that will constrain the range of possible temperature -pressure 
profiles and chemical abundances. The first model does not 
use a ramp model component and has the lowest value for the 
eclipse depth and the second model uses the sinusoidal ramp 
model that has the highest values for the eclipse depth. The 
SDNR and BIC values with the eclipse depths are given in 
Table 3. 

5.2. Channel 2-4.5 \im 

Channel 2 and 4 were observed at the same time. We 
first modeled each channel separately, determining the best 
aperture size, time-variability (ramp) model, and bin size for 
BLISS. Then we applied the joint fit. 




2.25 2.50 2.75 3.00 3.25 3.50 3.75 4.00 4.25 
Photometry Aperture Size (Pixels) 

Fig. 8. — Channel 2 comparison between linear and no ramp models. Plot 
shows SDNR and ABIC i'.v. aperture size. A lower value indicates a better 
model fit. 

The observation in channel 2 lasted 5.5 hours. There was 
no stabilization period observed in the data, so no points were 
removed from the analysis. The best aperture radius is 2.50 
pixels (see Figure 8). We also test the dependence of eclipse 
depth on aperture radius (see Figure 9). The eclipse depths 
are well within la. 
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Fig. 9. — The best-fit eclipse depths as a function of photomefi'y aperture 
size for wa014bs21. The red point indicates the best aperture size for that 
channel. The echpse-depth uncertainties are the result of 1000 MCMC itera- 
tions. The trend shows insignificant dependence of eclipse depth on aperture 
size (much less than Icr). 



Prior to the science observations in channels 2 and 4, we ob- 
served a 212-frame preflash, exposing the array to a relatively 
bright source to reduce the time-dependent sensitivity (ramp) 
effect in channel 4 (Charbonneau et al. 2005; Harrington et al. 
2007; Knutson et al. 2008). This produced saturation within 
the 42 minutes of the observation and led to a smaller rate of 
increase (see Figure 10). 
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Fig. 10. — Preflash light curve. The preflash source is bright compared to 
WASP-14b, which ramps airay sensitivity before the science observations. 
Without a preflash, similar observations generally show a steeper and longer 
ramp. 

We tested the full set of ramp equations and discarded ob- 
vious bad fits after shorter runs. Among acceptable ones, the 
lowest BIC value (see Table 4) determined that the final best 
ramp model is linear (see Eq. 5). 



R{t) = r,{t-Q.5)+l, 
where t is orbital phase and r\ is a free parameter. 



(5) 



TABLE 4 
Channel 2 Ramp Models 



Ramp Model 


SDNR 


BIC 


Eclipse 








Depth [%] 


No Ramp 


0.0044726 


2964.2 


0.224 ± 0.012 


Linear 


0.0044725 


2971.9 


0.224 ± 0.018 


Quadratic 


0.0044723 


2979.9 


0.241 ± 0.025 


Rising 


0.0044726 


2980.1 


0.224 ± 0.021 


Lin+Log 


0.0044690 


2983.9 


0.228 ± 0.017 



Each observation ended with a 10-frame, post-eclipse ob- 
servation of blank sky in the same array position as the sci- 
ence observations to check for warm pixels in the photometric 
aperture. There were none. 

To remove intrapixel variability we again apply our new 
BLISS technique. We also tested the intrapixel effect with a 
second-order, two-dimensional polynomial (see Eq. 4), con- 
firming that the BLISS technique gives a better fit (see Table 
2). Position precision in the channel 2 is 0.02 pixels for x and 
0.014 pixels for y. The best bin size is 0.028 pixels in x and 
0.023 pixels in y. The best aperture size, ramp model, and 
bin size for the intrapixel variability are then used in our joint 
fit (see Section 5.4), which gave us the eclipse depths and the 
brightness temperatures in Section 5.4. 

TABLE 5 
Channel 4 Ramp Models 



Ramp Model 


SDNR 


BIC 


Eclipse 








Depth [%] 


No Ramp 


0.0039799 


1459.2 


0.181 ± 0.013 


Linear 


0.0039770 


1464.3 


0.182 ± 0.012 


Rising 


0.0039799 


1466.4 


0.198 ±0.030 


Quadratic 


0.0039763 


1471.3 


0.181 ±0.018 


Lin+Log 


0.0039799 


1481.0 


0.181 ±0.024 



5.3. Channel 4 - 8.0 yun 

Again, no stabilization period was observed in this data set, 
hence no data points were removed from the 8.0 |^m dataset. 
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TABLE 6 

3-JoiNT Best-Fit Eclipse Light-Curve Parameters . 



Parameter 


Channel 1 


Channel 2 


Channel 4 


Array Position (x, pix) 


14.16 


23.82 


24.6 


Array Position (y, pix) 


15.69 


24.11 


21.9 


Position Consistency" {5y, pix) 


0.005 


0.02 


0.021 


Position Consistency" (Sy, pix) 


0.012 


0.014 


0.025 


Aperture Size (pix) 


2.75 


2.5 


3.5 


Sky Annulus Inner Radius (pix) 


8.0 


12.0 


12.0 


Sky Annulus Outer Radius (pix) 


20.0 


30.0 


30.0 


System Flux F, ( (jjy) 


102802.0 ±4.0 


66081.0 ±8.0 


24379.0 ± 4.0 


Eclipse Depth (%) 


0.187 ±0.007 


0.224 ±0.018 


0.1 82 ±0.022 


Brightness Temperature (K) 


2292.0 ±40.0 


2292.0 ±97.0 


1637.0 ± 120.0 


Eclipse Midpoint (orbits) 


0.4825 ± 0.0003 


0.4841 ±0.0005 


0.4841 ±0.0005 


Eclipse Midpoint (BJDuxc " 2,450,000) 


5274.6609 ± 0.0007 


4908.9288 ±0.0011 


4908.9288 ±0.0011 


Eclipse Midpoint (BJDxdb " 2,450,000) 


5274.6617 ±0.0007 


4908.9296 ±0.0011 


4908.9296 ±0.0011 


Eclipse Duration (f4_| , hrs) 


2.58 ±0.03 


2.58 ±0.03 


2.58 ±0.03 


Ingress/Egress Time (ti-i , hrs) 


0.291 ±0.007 


0.291 ±0.007 


0.291 ±0.007 


Ramp Name 


linear 


Hnear 


hnear 


Intrapixel Method 


BLISS 


BLISS 


- 


Ramp, Linear Term (ri ) 


0.0044 ±0.0010 


-0.0017 ±0.0031 


-0.005 ± 0.004 


Minimum Number of Points Per Bin 


4 


5 


- 


Total Frames 


13760 


2982 


1491 


Frames Used 


12308 


2944 


1433 


Rejected Frames (%) 


0.486919 


0.335345 


3.89001 


Free Parameters 


6 


4 


3 


Number of Data Points in Fit 


12308 


2944 


1433 


AIC Value 


12313.7 


4384.1 


4384.1 


BIC Value 


12358.2 


4428.79 


4428.79 


SDNR 


0.00331051 


0.0044727 


0.00397753 


Uncertainty Scaling Factor 


0.178797 


0.542751 


0.585052 


Photon-Limited S/N (%) 


72.7 


90.4 


68.1 



"RMS frame-to-frame position difference. 

Since we had the pre-calibration frames for this channel, 
which made the sensitivity increase more stable and less pro- 
nounced, the linear ramp was the best model for the observed 
time-variability (see Table 5). This is also confirmed with our 
best fit criteria and the lowest BIC value. No ramp model is 
excluded because channel 4 is known to exhibit strong ramp 
effect (Charbonneau et al. 2005; Harrington et al. 2007). 

The best aperture size is 3.50 pixels (see Figure 11). We 
again test for the dependence of eclipse depth with aperture 
size (see Figure 12). 





♦ ♦ No Ramp 

♦ ♦ Best aperture 











♦— ♦ Linear Ramp 
♦ ♦ Best aperture 










2.25 2.50 2.75 3.00 3.25 3.50 3.75 4.00 4.25 
Photometry Aperture Size (Pixels) 

Fig. 1 1 . — Channel 4 comparison between linear and no ramp models. The 
plot shows SDNR and ABIC vs. aperture size. A lower value indicates a 
better model fit. 

Even though intrapixel variability is not so strong in the 
Spitzer channels 3 and 4, channel 4 does exhibit a pixellation 
effect (see Stevenson et al. 2011 and Anderson et al. 2010), 
which justifies testing whether BLISS can give a better fit. 
Upon testing a full set of bin sizes, we concluded that NNI 



Photometry Aperture Size (Pixeis) 

Fig. 12. — The best-fit eclipse depths as a function of photometry aperture 
size for wa014bs41. The red point indicates the best aperture size for that 
channel. The eclipse-depth uncertainties are the result of 1000 MCMC itera- 
tions. This channel has the lowest S/N (~8). The aperture size of 2.25 pixels 
shows excess noise. Excluding it, the trend exhibits insignificant dependence 
of eclipse depth on aperture size (less than Icr). 



always outperforms BLI. That is the indication that intrapixel 
variability is very weak or non existing and can be ignored. 

The best aperture size and ramp model are used in our joint 
fit (see Section 5.4), which gave us the eclipse depths and the 
brightness temperatures in Section 5.4. 

5.4. Joint Fit 

To produce final joint-model fits, we run 10* MCMC itera- 
tions for all channels simultaneously, fitting a common eclipse 
duration for channel 1 , 2 and 4, and a common midpoint time 
for the channels observed together (channel 2 and 4). We 
again used the informative prior for ingress and egress times 
(Joshi et al. 2009 and David R. Anderson, personal commu- 
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TABLE? 
Best Joint-Model Fits 



Label 


Wavelength 


Ramp 


Intrapixel 


Bin Size 


Bin Size 


Eclipse 


Brightness 




[^m] 


Model 


Model 


in X 


in y 


Depth [%] 


Temp. [K] 


wa014bsll 


3.6 


lineal' (see also section 5) 


BLISS 


0.004 


0.01 


0.187 ±0.007 


2292 ± 40 


wa014bs21 


4.5 


linear 


BLISS 


0.028 


0.023 


0.224 ±0.018 


2292 ± 97 


wa014bs41 


8.0 


linear 








0.182 ±0.022 


1637 ± 120 




Fig. 13. — Left: Parameter correlations for wa014bsll. The background color depicts the absolute value of the correlation coefficient. Right: Parameter 
correlations for wa014bs21 and wa014bs41, produced in the joint fit. The background color depicts the absolute value of the correlation coefficient. Every lOO'*" 
step in the MCMC chain is used, to decoirelate consecutive values in both datasets. 
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Fig. 14. — Left: Parameter histograms for wa014bsll. Right: Parameter histograms for wa014bs21 and wa014bs4 1, produced in the joint fit. EverylOO step 
in the MCMC chain is used, to decorrelate consecutive values in both data sets. 



nication, 2010). The convergence was tested every 10^ steps. 
Our models converged after 5%, and the Gelman & Rubin 
(1992) diagnostic dropped below 1% for all free parameters. 
Figure 13 and Figure 14 show marginal posterior distributions 
(the parameter histograms) and plots of the pairwise correla- 
tions for channel 1 and joint channels 2 and 4. These confirm 
that the phase space minimum is global, and define the value 
of parameter uncertainties. The best-fit eclipse light curve pa- 
rameters from the joint fit and the best model fits are given in 



Table 6 and Table 7 respectively. 

6. ORBIT 

We fit the midpoint times from the Spitzer lightcurves si- 
multaneously with the available radial velocity curves and 
transit photometry, in order to provide updated measurements 
of system orbital parameters (see Table 8). The timing of sec- 
ondary eclipse is a strong constraint on the shape and orien- 
tation of the orbit. The two eclipses occur at phases 0.4825 
± 0.0003 and 0.4841 ± 0.0005 (using the Joshi et al. 2009 
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ephemeris), with a weighted mean after a 37 s Hght-time 
correction of 0.48273 ± 0.00025, indicating that ecoscj = - 
0.02713 ± 0.00039. The phases differ from each other by ap- 
proximately 3tT, but depend strongly on the accuracy of the 
ephemeris used to compute them. 



TABLE 8 
Eccentric Orbital Model. 



Parameter 


Value 




0.0831 ±0.0021 


ecos ix>^ 


-0.02557 ±0.00038 


e 


0.087 ±0.002 




-107.1 ±0.5 


P (days)" 


2.2437661 ±0.0000011 


To (MJD)" 


2454827.06666 ± 0.00024 




990 ± 3 


7 (ms"')" 


-4,987.9 ± 1.6 




162 


^Free parameter in MCMC fit 



We fit a Keplerian orbit model to our secondary eclipse 
times along with radial velocity data from Husnoo et al. 
(2011) and Joshi et al. (2009), and transit timing data from 
both amateur observers and WASP-14b's discovery paper 
(Joshi et al. 2009). The entire data set comprised 38 RV 
points, six of which were removed due to the Rossiter- 
McLaughlin effect, 30 transits, and two eclipses (see Table 
9). All times were adjusted to BJDjpg (Eastman et al. 2010), 
and 37 seconds were subtracted from the eclipse times to ap- 
proximately account for the eclipse-transit light-time. The 
errors were estimated using our MCMC routine. From this 
fit we determine that e = 0.087 ± 0.002 and w = 107.1 ± 
0.5°. With our new data, we refine the ephemeris to Tbjd-iob - 
2454827.06666(24)+2.2437661(ll)A', where T is the time of 
transit and is the number of orbits elapsed since the transit 
time (see Table 8). We find that the new ephemeris damp- 
ens the difference between the two eclipse phases to less than 
1.6(7. Performing an ephemeris fit to the transit and eclipse 
data separately shows that the transit and ecUpse periods dif- 
fer by (1.1 ± 0.8) X lO'^days, a 1.5(7 result that limits apsidal 
motion, oj, to less than 0.0024 °day~^ at the 3cr level (Gimenez 
& Bastero 1995). 

The results confirm an eccentric orbit for WASP- 14b and 
improve knowledge of other orbital parameters. 

7. ATMOSPHERE 

We modeled the dayside atmosphere of WASP- 14b using 
the exoplanetary atmospheric modeling method developed by 
Madhusudhan & Seager (2009, 2010). We use a ID hne-by- 
line radiative transfer code to model the planetary atmosphere 
under the assumption of local thermodynamic equiUbrium 
(LTE), hydrostatic equilibrium, and global energy balance at 
the top of the atmosphere. The latter condition assumes that 
the integrated emergent planetary flux balances the integrated 
incident stellar flux, accounting for the Bond albedo (Ab) and 
possible redistribution of energy onto the night side. Our 
model uses parametrized prescriptions to retrieve the temper- 
ature structure and chemical composition from the observa- 
tions, as opposed to assuming radiative and chemical equilib- 
rium with fixed elemental abundances (Burrows et al. 2008; 
Fortney et al. 2008). 

We explore the model parameter space in search of the 
best-fitting models for a given data set. The model param- 
eterization is described by Madhusudhan & Seager (2009, 



TABLE 9 
Transit Timing Data. 



Mid-Transit Time (BJDxdb) 


Uncertainty 


Source" 


2455695.4082 


0.0012 


V. Slesarenkno, E. Sokov ^ 


2455668.4790 


0.0011 


Frantisek Lomoz 


2455652.7744 


0.0014 


Stan Shadick, C. Shiels'' 


2455650.5307 


0.0018 


Lubos Brat 


2455650.52789 


0.00076 


Martin VraSt'Sk 


2455650.52566 


0.00067 


Jaroslav Tmka 


2455632.5807 


0.0011 


E. Sokov, K. N. Naumov ^ 


2455318.45101 


0.00085 


Anthony Ayiomamitis 


2455302.7464 


0.0010 


Stan Shadick^ 


2455264.6021 


0.0012 


Hana Kucakova 


2455264.6017 


0.0013 


Radek Kocian*^ 


2455219.7290 


0.0012 


Lubos Brat 


2454979.643 


0.003 


Wiggins, AXA 


2454968.426 


0.001 


Srdoc, AXA 


2454950.4831 


0.0021 


Jesionldewicz, AXA 


2454950.4746 


0.0014 


Lubos Brdt 


2454950.4745 


0.0018 


Hana Kucdkovi 


2454950.4731 


0.0021 


Pavel Marek 


2454950.4728 


0.0014 


Wardak, AXA 


2454943.7427 


0.0006 


Dvorak, AXA 


2454941.49799 


0.00081 


Jaroslav Tmka 


2454941.4916 


0.0019 


Frantisek Lomoz 


2454934.765 


0.001 


Brucy Gary, AXA 


2454932.5246 


0.0014 


Radek Df eveny 


2454932.5232 


0.0011 


Lubos Brit 


2454932.5222 


0.0013 


Jaroslav Tmka 


2454932.5219 


0.0015 


T. Hynek, K. Onderkovd 


2454914.5753 


0.0008 


Naves, AXA 


2454887.6457 


0.0014 


Georgio, AXA 


"The Amateur 


Exoplanet 


Archive (AXA), 



http://brucegary.net/AXA/x.htm) and Transiting ExoplanetS and Can- 
didates group (TRESCA), http://var2.astro.cz/EN/tresca/index.php) 
supply their data to the Exoplanet Transit Database (ETD), 
http://var2.astro.cz/ETD/) which performs the uniform transit anal- 
ysis described by Poddany et al. (2010). The ETD web site provided 
the AXA and TRESCA numbers in this table, which were converted to 
BJDtdb- 

''Physics and Engineering Physics Department, University of 
Saskatchewan, Saskatoon, Saskatchewan, Canada, S7N 5E2 

'^KociSn R., Johann PaUsa, Observatory and Planetarium, Technical 
University Ostrava, 17. Ustopadu 15, CZ-708 33 Ostrava, Czech Republic 

''Sokov E., Naumov K., Slesarenko V. et al., Pulkovo Observatory of 
RAS, Saint-Petersburg, Russia 

■^Municipal Observatory in Slany Czech Republic 

^Project Eridanus, Observatory and Planetarium of Johann Palisa in Os- 
trava 

2010). The sources of opacity in the model include molec- 
ular absorption due to HjO, CO, CH4, COj, TiO, and VO, 
and collision-induced absorption (CIA) due to H2-H2. Our 
molecular line lists are obtained from Freedman et al. (2008), 
Freedman (personal communication, 2009), Rothman et al. 
(2005); Karkoschka & Tomasko (2010), and Karkoschka (per- 
sonal communication, 2011). Our CIA opacities are obtained 
from Borysow et al. (1997) and Borysow (2002). We explore 
the model parameter space using a MCMC scheme, as de- 
scribed by Madhusudhan & Seager (2010). However, since 
the number of model parameters (A^ = 10) exceed the number 
of data points {Ndata = 3), our goal is not to find a unique fit 
to the data but, primarily, to rule out regions of model phase 
space that are excluded by the data. In order to compute the 
model planet-star flux ratios to match with the data, we di- 
vide the planetary spectrum by a Kurucz model of the stellar 
spectrum derived from Castelli & Kurucz (2004). Our mod- 
els allow constraints on the temperature structure, molecular 
mixing ratios, and a joint constraint on the albedo and day- 



Thermal Emission of WASP- 14b 



11 



WASP- 14b 




10 






IRAC 8.0 


\ c 


IRAC 5,8 ■ 




IRAC 4.5 


'\ 


IRAC 3.6 










V. ■- 













1.0 1.5 2.0 2.5 J.D 0.2 0.4 0.5 0.8 1.0 
Normalized Contribution 



T (1000 K) 



0.4 0.6 o.a 1 
nolized Contribution 



Fig. 15. — Left: Observations and model spectra for dayside emission from WASP-14b. The blue filled circles with error bars show our observations in Spitzer 
channel 1 (3.6 |J.m), 2 (4.5 ^m), and 4 (8.0 [im). For the 3.6 |a.m channel, two values are shown, in blue and brown, corresponding to different ramp models 
used in deriving the eclipse depths (see Section 5.1). The green and red curves in the main panel show two model spectra without thermal inversions and with 
different chemical composition; the corresponding pressure-temperature (P—T) profiles are shown in the inset. The green and red circles show the model spectra 
integrated in the Spitzer IRAC bandpasses. Both profiles lead to good fit to the observations, but with different molecular compositions. The green model has 
molecular abundances in chemical equilibrium assuming solar elemental abundances. The red model has 10 times lower CO and 6 times higher H2O compared to 
solar abundance chemistry, and provides a better fit to the data compared to the green model. The black dotted lines show three blackbody planet spectra at 1500 
K, 2200 K, and 2600 K. Right: P-T profiles and contribution functions corresponding to the models shown on the left panel. Panel a shows two P—T profiles 
that correspond to the P-T profiles in the inset figure, left panel. Solid line corresponds to the green line, and dotted line to the red line. Panel b and c show the 
contribution functions in the four Spitzer channels for each model. The legend shows the channel center wavelength in ^.m and the curves are color-coded by the 
channel. All the contribution functions are normalized to unity. 



night redistribution. 

We find that strong constraints can be placed on the pres- 
ence of a thermal inversion in WASP- 14b even with our cur- 
rent small set of observations. At an irradiation of 3 x lO"^ 
ergss 'cm^^, WASP- 14b falls in the class of extremely irradi- 
ated planets that are predicted to host thermal inversions ac- 
cording to the TiOA'O hypothesis of Fortney et al. (2008). 
However, the present observations do not show any distinct 
evidence of a thermal inversion in the dayside atmosphere of 
WASP-14b. While the brightness temperatures in the 3.6 and 
4.5 [im channels are consistent with a blackbody spectrum of 
the planet at T^2200 K, the 8 micron flux deviates substan- 
tially from the assumption of a blackbody with a brightness 
temperature of 1637 ± 120 K. In the presence of a thermal 
inversion, the flux in the 8 |J.m channel is expected to be much 
higher than the fluxes in the 3.6 and 4.5 |J.m channels due 
to emission features of water vapor, and methane, if present. 
The low flux observed at 8 |a.m, therefore, implies strong wa- 
ter vapor and/or methane in absorption, implying the lack of 
a significant temperature inversion (see Madhusudhan & Sea- 
ger 2010 for a discussion on inferring thermal inversions). We 
also note that, as mentioned in Section 5.1, the observations 
yield different planet-star flux contrasts in the 3.6 [im channel 
for different choices of ramp models. However, as shown in 
Figure 15, the two extreme values are still consistent at the 
l-cr level, and as such, lead to similar model conclusions. 

The observations provide only marginal constraints on the 
chemistry, although there are hints of enhanced H2O over CO. 
We find that the observations can all be explained by compo- 
sitions within a factor of 10 of those obtained from thermo- 
chemical equilibrium assuming solar abundances (TEj^i^^). 
Figure 15, shows two model spectra along with the observa- 
tions, one with molecular abundances at TEj.(,,2j., and another 
with 10 times lower CO and 6 times higher H2O. The latter 
model fits the data marginally better. A slightly lower CO 
is favored because of the slightly higher 4.5 \im flux com- 
pared to the 3.6 |J.m flux, which means lower absorption due 
to CO. Higher absorption due to H2O is favored by the low 8 



|j.m point. In principle, a lower CO and a higher H2O, com- 
pared to TEjqi^j values are both possible by having a lower 
C/O ratio compared to the solar C/O ratio (0.54). However, 
more data would be required to confirm the low CO require- 
ment, because a blackbody of ^2200 K fits the 3.6 and 4.5 
|j.m points just as well. We have also considered models with 
high C/O ratios (e.g. C/O > 1), which can lead to strong CH4 
absorption in the 3.6 |J.m and 8 |J.m channels (e.g., Madhusud- 
han et al. 201 la,b), instead of H2O absorption in the low C/O 
models shown in Figure 15. However, we find that the amount 
of methane possible in high C/O models of WASP- 14b is in- 
adequate to cause the strong absorption in the 3.6 [im and 8 
|j.m channels required by the data; note that in high C/O mod- 
els H2O is depleted and CH4 is enhanced. This would be a 
good motivation for future observations in the near IR, from 
ground and space, to place further constraints on the temper- 
ature structure and composition, especially the C/O ratio, of 
the dayside atmosphere of WASP- 14b. 

The models explaining the observations require relatively 
low day-night energy redistribution in WASP- 14b. The large 
brightness temperature in the 3.6 [im channel, which probes 
the deep atmosphere at pressures of ^1 bar, is indicative of a 
hot planet photosphere. For zero-albedo, we find that about 
^0.25 of the energy incident on the dayside can be redis- 
tributed to the night side. For non-zero albedos this fraction is 
even lower; since the quantity we constrain is 77 = (1-Ab)(1- 
fj), where Ag is the Bond albedo, and f^. is the fraction of the 
dayside incident energy redistributed to the night-side. Our 
finding of low redistribution can be verified by thermal light 
curves in the future (e.g., Knutson et al. 2007, 2009a). 

8. DISCUSSION 

The absence of a thermal inversion in the dayside at- 
mosphere of WASP- 14b provides important constraints on 
inversion-causing phenomena in irradiated atmospheres. The 
canonical argument for thermal inversions in hot Jupiter at- 
mospheres is via TiO and VO. Hubeny et al. (2003) and Fort- 
ney et al. (2008) suggested the possibility that TiO and VO 
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could exist in gas phase in high-temperature atmospheres - 
the TiO/VO hypothesis. Being strong absorbers of incident 
visible Ught, TiO and VO could cause thermal inversions in 
hot Jupiter atmospheres. On the other hand, Spiegel et al. 

(2009) showed that the high mean molecular masses of TiO 
and VO would lead to significant gravitational settling of 
these molecules, thereby depleting them from the upper at- 
mospheres, unless strong vertical mixing keeps them aloft. 
Additionally, the abundances of inversion-causing molecules 
might also be influenced by stellar activity and photochem- 
istry (Knutson et al. 2010). Consequently, the real cause of 
thermal inversions in irradiated atmospheres is currently un- 
known. Nevertheless, models used to infer thermal inversions 
in the literature have either used parametrized visible opacity 
sources (Burrows et al. 2008) or parametric temperature pro- 
files (Madhusudhan & Seager 2009, also used in the present 
work). To first order, the lack of a thermal inversion in WASP- 
14b might indicate that the vertical mixing in the day side at- 
mosphere of WASP- 14b is weaker compared to the downward 
diffusion of TiO and VO. 

Until now, Spitzer observations have provided us with a 
range of hot Jupiter planets that are strongly irradiated by 
their parent stars, with brightness temperatures between 1000- 
2000 K. The inferences of thermal inversions from thermal 
spectra of hot Jupiters result from observing flux excesses in 
molecular bands where strong absorption is expected (Mad- 
husudhan & Seager 2010). In principle, detection of a ther- 
mal inversion in high S/N observations is possible with just 
two Warm Spitzer channels if there is a high difference in 
fluxes in channels 1 and 2. Similar conclusions have been 
reached based on the flux excesses in certain Spitzer channels 
by e.g., Knutson et al. (2008, 2009a); Madhusudhan & Seager 

(2010) ; Machalek et al. (2009); O'Donovan et al. (2010) and 
Christiansen et al. (2010). Based on the TiO/VO hypothesis 
described above, Fortney et al. (2008) suggested that depend- 
ing on the level of irradiation from their parent star, irradiated 
planets can fall into two categories: the very highly irradiated 
atmospheres that host thermal inversions and the less irradi- 
ated ones that do not. However, recent observations have re- 
vealed several counter-examples to this hypothesis. Machalek 
et al. (2008, 2009) present evidence for a temperature inver- 
sion in XO-lb, despite low irradiation of the planet (Teq = 
1209 K), while Fressin et al. (2010) show no thermal inver- 
sion, although TrES-3b is a highly irradiated planet (Teq = 
1643 K). Similarly, WASP-12b, one of the most irradiated hot 
Jupiters known, has also been reported to lack a significant 
thermal inversion (Madhusudhan et al. 2011a). In this paper, 
we present WASP- 14b as an another example that demon- 
strates that planets studied to date cover a wide range and do 
not fall distinctively into two categories (Gillon et al. 2010). It 
is possible, that additional parameters (e.g., metalUcity or sur- 
face gravity) influence the presence or the absence of a tem- 
perature inversion. However, more observations are needed to 
explain the lack of the inversion in this planet. 

9. CONCLUSIONS 

During two secondary ecUpse events, Spitzer observed 
WASP-14b in the three IRAC channels: 3.6, 4.5, and 8.0 |xm. 
All eclipses have a high S/N ratio (3.6 [im channel ^25, 4.5 
|xm channel ^12, 8.0 ^im channel ^8), which allowed us to 
constrain the planetary spectrum and orbital parameters. 

Our observations in the three IRAC channels probe the at- 
mosphere at pressures between 0.01 and 1 bar and indicate the 



absence of a significant thermal inversion in the dayside atmo- 
sphere of WASP- 14b. The lack of a thermal inversion in the 
highly irradiated atmosphere of WASP- 14b contradicts pre- 
vious predictions that the most irradiated hot Jupiters should 
have thermal inversions due to gaseous TiOA^O (Fortney et al. 
2008). Additionally, our observations place nominal con- 
straints on the chemical composition and day-night energy 
redistribution in the atmosphere of WASP- 14b. We find that 
the data can be explained by non-inversion models in chemi- 
cal equiUbrium with nearly solar abundances. A factor of 10 
less CO and a factor of 6 higher HjO, compared to those ob- 
tained with solar abundances, explain the data to within the 
1(7 uncertainties, on average. Such CO depletion and H2O en- 
hancement are, in principle, possible in chemical equilibrium 
with C/O ratios lower than solar. More data are required to 
better constrain the atmospheric composition of WASP- 14b. 

Because the planet is much brighter than its predicted equi- 
librium temperature for uniform redistribution (Teq = 1 866^"), 
the best-fitting models require relatively low day-night energy 
redistribution in WASP- 14b. We find an upper limit on the 
fraction of energy redistributed to the night side to be 0.25 
(for zero Bond albedo). 

WASP- 14b is one of the most massive transiting planets 
known along with CoRoT-3b (Triaud et al. 2009; Deleuil et al. 
2008), HAT-P-2b (Bakos et al. 2007; Winn et al. 2007; Loeil- 
let et al. 2008b), XO-3b (Hebrard et al. 2008; lohns-Krull 
et al. 2008; Winn et al. 2008), and WASP-18b (Nymeyer et al. 
2011). With the exception of WASP-18b, all of these objects 
have large eccentric orbits. Classic theories suggest that there 
is a trend for the inner planets to have circular orbits due to 
tidal orbital decay. At distances a < 0. 1 AU any orbital eccen- 
tricity should induce strong tidal energy dissipation, leading 
to circular orbits in relatively short time scales (typically a few 
Myr, compared to common system ages of few Gyr). How- 
ever, Pont et al. (2011) argue that the circularization time- 
scale scales with the planet-star mass ratio, and is also a steep 
function of the orbital separation, scaled to the planet radius 
(see their Figure 3). For planets with M > Mj, the mass-period 
relation (see their Figure 2) suggests that heavier planets get 
circularized very close to their parent star, or they do not ever 
reach circularization in their lifetime. A possible explanation 
is that the planet raises tides on its host star strong enough that 
the angular momentum of the planet is transfered to the stel- 
lar spin, and the planet gets swallowed by the star. This does 
not oppose the classical tide theory (e.g., Goldreich & Soter 
1966), rather suggests that stopping mechanisms and tidal cir- 
cularization are related. WASP-14b also has unusually high 
density for hot Jupiters, similar to that of some rocky planets 
(4.6 gcm"^). Having a strong signal makes this planet ideal 
for further exploration of this topic and of possible formation 
mechanisms of similar objects. 
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Thermal Emission of WASP- 14b 
A. SYSTEM PARAMETERS 



TABLE 10 
System Parameters of WASP- 14. 



PARAMETER 


VALUE 


REFERENCE 


Eclipse Parameters 


Eclipse JMidpoint (BJDxdb) (2010-03-18) 2455274.6617 ± 0.0007 


a 


EcUpse Midpoint (BJDtdb) (2009-03-18) 2454908.9296 ± 0.001 1 


a 


Eclipse Duration f4_i (hrs) 


2.58 it 0.03 


a 


Transit depth, {RplRstar)^ 


0102 ^o-^^ 

u.uiuz _Q 0003 


c 


Ingress/Egress Time f2-i (hrs) 


0.291 ± 0.007 


a 


Orbital Parameters 


Orbital Period, P (days) 


2.2437661 ± 0.0000011 


b 


Semimajor axis, a (AU) 


0.036 ± 0.001 


c 


Transit Time (BTD^qb) 


2454827.06666 ± 0.00024 


b 


Orbital eccentricity, e 


0.087 ± 0.002 


b 


Argument of pericenter, ix> (deg) 


-107.1 ± 0.5 


b 


Velocity semiamplitude, K (ms"') 


990.0 ± 3 


b 


Centre-of-mass velocity 7 (ms"') 


-4987.9 ± 1.6 


b 


Stellar Parameters 


Spectral Type 


F5V 


c 


Mass, Mt (Mq) 
Radius, R, (Rq) 
Mean density, p* (pq) 


1 211 +"'27 
-0.122 

1 306 

I.OUD _Q 073 
542 


c 
c 
c 


Effective temperature, Tgff (K) 
Surface gravity, log g* (cgs) 


6475 ± 100 

A OOT +0.043 
H..Z.O/ _o.o38 


c 
c 


Projected rotation rate, v« sin(i) (kms" ) 4.9 ±1.0 


c 


Metallicity [M/H] (dex) 


0.0 ± 0.2 


c 


Age (Gyr) 


~0.5-1.0 


c 


Distance (pc) 


160 ± 20 


c 


Lithium abundance, log N(Li) 


2.84 ± 0.05 


c 


Planetary Parameters 


Mass, Mp (Mj) 


7 041 +0.508 


c 


Radius, R,, (Rj) 


1 901 +0.075 
-0.082 


c 


Surface gravity, log gp (cgs) 


4 010 

t.UlU _o 042 


c 


Mean density, pp (gem ) 


4.6 


c 


Equilibrium temperature (A=0), T, 


,,(K) 1866.12 


c 



a - our analyses (parameters derived using 3-joint fit, see Section 5.4), 
b - our analyses (see Section 6) 
c - Joshi et al. (2009), 



